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We show that the magnetoplasmon collective modes in quasi-two-dimensional semiconductor 
quantum dots can be parametrically amplified by periodically modulating the magnetic field per- 
pendicular to the nanostructure. The two magnetoplasmon modes are excited and amplified simul- 
taneously, leading to an exponential growth of the number of bosonic excitations in the system. We 
further demonstrate that damping mechanisms as well as anharmonicities in the confinement of the 
quantum dot lead to a saturation of the parametric amplification. This work constitutes a first step 
towards parametric amplification of collective modes in many-body fermionic systems beyond one 
dimension. 
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I. INTRODUCTION 

The swing is the paradigm for the phenomenon of para- 
metric resonance: 1 The periodic motion of the swinger's 
legs leads to a periodic modulation of the effective length 
of the pendulum, hence of its frequency. This results in 
the exponential amplification of the motion if the period 
of the modulation is chosen to be commensurate with the 
natural frequency of the pendulum. Quantum mechani- 
cally, if one considers a single harmonic oscillator whose 
frequency is periodically varied in time, the parametric 
modulation of the Hamiltonian leads to an exponential 
divergence of the corresponding lowering and raising op- 
erators, hence increasing the number of bosonic quanta 
in the system. As in the classical case, 1 this effect is most 
efficient if the pumping frequency is twice as much as the 
natural frequency of the oscillator. For a parametrically 
modulated many-body bosonic system, only the eigen- 
modes fulfilling the resonance condition are amplified. 2 ' 3 
Parametric resonance is thus an interesting and valuable 
spectroscopic tool for many-body quantum systems, es- 
pecially in the context of cold atomic gases 4 where the 
tunability of the experimental setups enables one to para- 
metrically amplify Bosc-Einstein condensates. 5,6 

In this paper, we address the question of parametric 
resonances in many-body fermionic systems. In con- 
trast to the bosonic case, any mean-field treatment of 
the interactions between fermions preserves the fermionic 
nature of the quasiparticles: The parametric amplifi- 
cation of fermionic modes is thus blocked due to the 
Pauli principle. However, many-body fermionic systems 
are known to exhibit bosonic collective modes such as, 
e.g., plasmons (charge density waves) and magnons (spin 
waves). 7 ' 8 It has been shown in various contexts that 
such collective modes can be parametrically excited in 
one-dimensional correlated fermionic systems. 9-12 Using 
the Luttinger liquid picture and bosonization techniques 
which enable one to treat the interactions among the 
particles (almost) exactly, 8 ' 13 it was proposed that the 
parametric modulation of the optical lattice where cold 



fermionic atoms are trapped could serve as a probe for 
the well-known spin-charge separation inherent to onc- 
dimcnsional correlated fermionic systems. 10,11 

Extending these ideas to two- and three-dimensional 
many-body fermionic systems is a formidable theoreti- 
cal challenge, as one needs to capture the relevant cor- 
relations by treating interactions beyond mean field. In 
this paper, we focus on the case of confined many-body 
fermionic systems, where collective modes exist due to 
the finite size of the system. 14 Examples of such collec- 
tive modes are plasmon excitations which correspond to 
the motion of the electronic center of mass in metallic 
nanoparticles 14-16 and in quasi-two-dimensional semicon- 
ductor quantum dots. 17-19 

It is well known that a static magnetic field 
induces field-dependent collective modes termed 
magnetoplasmons. 1 ' 18 These modes appear in partic- 
ular in quasi-two-dimensional semiconductor quantum 
dots, 20-22 and were recently investigated theoretically in 
metallic nanoparticles. 16 The magnetoplasmons, whose 
classical motion corresponds to an excitation of the 
electronic center of mass perpendicular to the magnetic 
field rotating clockwise and counterclockwise (see the 
inset in Fig. 1), have distinct frequencies that depend 
on the applied static magnetic field (see Fig. 1). In this 
work, we propose to trigger the parametric amplification 
of the magnetoplasmons by periodically modulating 
the magnetic field perpendicular to the nanostructure. 
Specifically, we focus on quasi-two-dimensional semi- 
conductor quantum dots. We show that the parametric 
modulation of the magnetic field leads to an exponential 
growth of the number of bosonic modes in the system. 
Wc further demonstrate how damping mechanisms and 
anharmonicities in the confinement of the quantum dot 
lead to the saturation of the parametric amplification. 

This article is organized as follows: We start in Sec. 
II by presenting our model of a two-dimensional quan- 
tum dot subject to a perpendicular magnetic field which 
is periodically modulated in time. In Sec. Ill, we show 
that such a modulation leads to the parametric ampli- 
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FIG. 1. (Color online) Frequencies oj+ and w_ [cf. Eq. (14)] of 
the two magnetoplasmon modes (rescaled by the confinement 
frequency ujo) as a function of the cyclotron frequency uj c = 
eBo/m. The pumping frequency exactly at resonance, Q = 
OJ+ + uj-, is also shown as a dotted line. Inset: Sketch of 
the two independent center-of-mass motions associated to the 
surface magnetoplasmon modes perpendicular to the static 
magnetic field B . 



fication of the magnetoplasmon collective modes. The 
effect of damping and anharmonicities on the parametric 
resonance of the magnetoplasmons is studied in detail in 
Sec. IV. We discuss in Sec. V the experimental conse- 
quences of our proposal before we conclude in Sec. VI. 
Several technical aspects of our work are presented in the 
appendices. 



II. MODEL 

We consider 7V C interacting electrons with effective 
mass m* and charge — e confined in a two-dimensional 
quantum dot. The quantum dot is subject to a spa- 
tially homogeneous perpendicular magnetic field B(<) = 
B(t)e z with 



B{t) = B [l +7? sin (fit)] 



which is periodically modulated at the frequency 57. In 
what follows, we assume the strength of the periodic 
modulation -q to be much smaller than one. The Hamil- 
tonian of the system reads 



ff (*)=E i(pi+ eA ( r - i )) 2+c/ ( r 

i=l L 

E^- 



(2) 



i,j=X 



with r s ; = Xie x + i/i&y the position of the ith electron, 
Pi its momentum and A(r, t) the time-dependent vector 



potential. Notice that we do not consider the spin de- 
gree of freedom as it is not relevant for understanding 
the parametric amplification of the magnetoplasmons, as 
the latter only involves the orbital degrees of freedom. 
Hence, we neglect the Zeeman term (which, for confined 
semiconductor quantum dots is relatively weak due to 
the small value of the associated g factor) and the spin- 
orbit coupling in Eq. (2) . The single-particle confinement 
U (r) appearing in Eq. (2) is approximated by a parabolic 
potential with confining frequency luq. 17,18 Including also 
quartic corrections that will be relevant for the following 
analysis (see in particular Sec. IV), we write it as 



U(r) 



04 4 
V 

4 ' 



(3) 



with ci4 > 0. Finally, V c . c in Eq. (2) stands for the 
electron-electron interaction. In the symmetric gauge 



. , x Bit) , 
A ( r , 1) = -—{-ye x + xe y ), 



(4) 



Pi , TTf- \ , m*LU 2 (t) 



one can rewrite the Hamiltonian (2) as 

tf(*)=E 

i 



r,l), 



(5) 



with l z _i = x Pi) z the z-component of the angular 
momentum, us(t) = eB(t)/m* = uj c [1 + rj sin (fit)}, and 
lu c = eBo/m* the cyclotron frequency. As we will show 
in Sec. Ill, the diamagnetic term oc w 2 (t) in Eq. (5) is 
responsible for the parametric amplification of the mag- 
netoplasmon collective modes. 

Introducing the collective coordinate for the electronic 
center of mass R = Xe x + Ye y = '^2 i Yi/N and its 
conjugated momentum P = Px&x + Py^y = S»P»j as 
well as the relative degrees of freedom r[ = — R and 
p' = pi — P/iV, 15 ' 16 the Hamiltonian (5) separates into 



H(t) = H^t) + H iel (t) + H t 

with the center-of-mass Hamiltonian 
P 2 M ( , u 2 {t) 



(1) H cm {t) = — + - 



2M ' 2 V" ° ' 4 J" ' 2 

(7) 

where M = N c m* , A 4 = N c a A , and L z = XP Y - YP X - 
The Hamiltonian for the relative coordinates reads 



If 



L 2 



(6) 



T ' 



rel 



(*)=E 



'2 



P: 

2 m 



2 z < 



m*uj 2 (t) 



,'2 



(8) 



Finally, the quartic part of the single-particle confine- 
ment (3) leads to the coupling between the center of mass 
and the relative coordinates in Eq. (6), with Hamiltonian 



H r = 



a 4 \ -> 



(|R- 



-r'| 4 -i? 4 -r' 4 



(9) 
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Notice that if one assumes the single-particle confine- 
ment (3) to be purely harmonic (04 = 0), there is no 
coupling in Eq. (6) between the electronic center of mass 
and the relative coordinates, as requested by the general- 
ized Kohn's theorem. 17,18 ' 23,24 The coupling (9) induces 
the damping of the center-of-mass excitations, and leads 
to a finite linewidth of the magnetoplasmon lines in far 
infrared absorption spectroscopy experiments. 17,18 

The time-independent, quadratic part of the center- 
of-mass Hamiltonian (7) can be diagonalized by means 
of Fock-Darwin states in terms of two magnetoplasmon 
excitations. 17 Introducing the variable £ = (X + \Y)/\/2 
and the bosonic operators associated to the magnetoplas- 
mons 



6+ = 



1 

71 



1 (a 



V2 v 



(10a) 



-£ 



with the oscillator length 



MK 2 + Wc 2 /4) 1/2 ' 



(10b) 



(11) 



we rewrite Eq. (7) as 

H cm (t) = [1 + <^/wi ( btb„ 



<j=± 

huj, 



+ n ^6f(t)(b\bt + b-b + ). (12) 



In Eq. (12), 



V ujI + w c 2 /4 



(13) 



and we introduced the dimensionless parameter a = 
Ait A jhbjQ <C 1. The frequencies of the two magneto- 
plasmon collective modes read 17 



4 



± 



(14) 



and are presented in Fig. 1 as a function of the cyclotron 
frequency w c . Without parametric modulation (77 = 0), 
the center-of-mass motion associated with the magneto- 
plasmon with frequency lu + rotates counterclock- 
wise (clockwise) in the plane perpendicular to the mag- 
netic field (see the inset in Fig. 1). Notice that one 
has uj+ > uj- due to the fact that the Lorentz force 
— iV e eR x Bo increases (decreases) the effective harmonic 
confinement seen by the collective mode with frequency 



The last symmetry breaking term cx +b + b_ in the 
Hamiltonian (12) is triggered by the parametric modula- 
tion of the magnetic field [see Eq. (13)]. As we will show 
in Sec. Ill, this term is responsible for the parametric am- 
plification of the two magnetoplasmon collective modes, 
provided that the pumping frequency is close to the 
resonance condition Q = uj + + cj_ (sec dotted line in Fig. 
1). The quartic corrections proportional to the parame- 
ter a in the Hamiltonian (12) represent a residual inter- 
action between the two bosonic modes, and lead to the 
saturation of the parametric amplification (see Sec. IV). 

In the absence of quartic corrections (a = 0), the 
Hamiltonian (12) is similar to the one encountered in 
the quantum theory of parametric amplifiers in quan- 
tum optics. 25-29 In this context, a pump laser interacts 
with a nonlinear crystal. Due to the second-order sus- 
ceptibility of the nonlinear media, a photon of the pump 
laser at frequency 17 splits into two photons at frequen- 
cies u + and cj_, a phenomenon called parametric down 
conversion. In a parametric amplifier, the two signals 
at frequencies lu+ and w_ are amplified by pumping the 
crystal at fl = uj + + w_ . 



III. PARAMETRIC RESONANCE OF THE 
MAGNETOPLASMONS 

We start our study of the possibility of parametrically 
amplifying the magnetoplasmon collective modes by ne- 
glecting the quartic part of the center-of-mass Hamilto- 
nian, setting a = in Eq. (12). Furthermore, we neglect 
all other sources of dissipation on top of the one induced 
by the quartic confinement. (For the simpler case of a sin- 
gle harmonic oscillator with periodically-modulated fre- 
quency, see Appendix A.) The Heisenberg equations of 
motion for the bosonic operators appearing in the center- 
of-mass Hamiltonian (12) thus read 



b± -- 
Introducing 



-iw± [1 ± Sf(t)] b ± 



^Sf(t)b 



(15) 



b ± (t) = exp I -i« ± ^ ds [1 ± 5f(s)}\ b ± (t), (16) 

Eq. (15) transforms, to first order in the small parameter 
77, into 

.Wc 

y 

The main contribution to the time evolution of these op- 
erators, within the rotating wave approximation, comes 
from the terms close to resonance, i.e., ~ ui+ +cj_ [cf. 
Eq. (13)], yielding 



b ± = -^Sf{t)e^+^b%. 



(17) 



b± ~ ie-'C 1 -^-"-)'^ 

where we introduced the "driving" frequency 

(wc/2) 2 



e = T) 



Vw 2 +w c 2 /4' 



(18) 



(19) 
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Equation (18) can be easily decoupled, yielding 



and 



b± + i(0 - lu + - uj-)b ± - — b± = 0. (20) 

Using the initial conditions b±(0) = 6±(0) and b±(0) = 
eb^-(0)/2 [cf. Eq. (18)], we finally obtain the solutions 



6±(t) 



(0_e i: 



a , t 



+ e in - 4 ) b ± {0) 



i £( e ^+*-e in - t )6^(0) 



with 



(21) 



(22) 



We thus see that tuning the pumping frequency O in the 
interval |0 — cj + — u—\ < e, the frequencies 0± acquire 
an imaginary part, leading to the parametric amplifica- 
tion of the bosonic magnetoplasmon collective modes. 
(For a classical treatment of the parametric resonance 
of the magnetoplasmons and specifically of the associ- 
ated classical trajectories of the electronic center of mass, 
see Appendix B.) Indeed, close to resonance, i.e., for 
|0 — uj+ — cj_| <C e, Eq. (21) simplifies to 



b±(t) 



= e~ iui±t 



cosh ( — ) b± (0) + sinh 



bUo) 



(23) 

showing that the two bosonic modes are exponentially 
amplified at a rate given by the frequency e of Eq. (19). 
On the contrary, off resonance (10 — uj+ — cj_| e), 
the bosonic operators evolve in time according to their 
unitary evolution, 



b±(t) 



^b±(0). 



(24) 



Notice that the exponential amplification of the operators 
(21) is a direct consequence of their bosonic nature. If the 
operators b± were fermionic, the Pauli principle would 
prevent the parametric amplification from occurring (see 
Appendix C for details). 

The center-of-mass subsystem is, prior to the para- 
metric modulation that starts at t = 0, in thermal 
equilibrium at the temperature T and thus described 
by the density matrix po — e~ Hcm (°^ kBT jZ with Z = 
tr{e~ ffcm (°)/ fcBT } the canonical partition function. Intro- 
ducing the average number of bosons in each mode 



N ± (t) = (bib ± )(t), 

as well as the "anomalous" average 



(25) 



<p(t) = (b+b-)(t) (26) 
where (0)(t) = tv{p O(t)}, we obtain using Eq. (23) 



N±(t) = ^{ [1 + n B {uj + ) +n B (w_)]cosh(ei) 
+ n B (w±) - 1 - n B (uj T )} 



(27) 



¥>(*) = 



,-iOt 



[l + n B (w+) + nB(w-)]sinh.(rf) (28) 



when the pumping frequency is close to resonance. In 
Eqs. (27) and (28), we introduced the Bose distribution 
n B (w) = [exp (hw/k-gT) — 1] _1 . Due to the fact that the 
magnetoplasmons are purely harmonic excitations of the 
electronic center of mass, (R)(t) — 0, the fluctuations 
of the center-of-mass coordinate (5R 2 )(t) = (R 2 )(t) — 
(R) 2 (t) reduce to 

(SR 2 )(t) = £ 2 [N+(t) + N_(t) + 1 + <p(t) + <p*(t)}. (29) 

With Eqs. (27) and (28), we obtain 

(SR 2 )(t) =l 2 [l + n B (u + ) + n B (w-)] 



x [cosh (et) + cos (fit) sinh (et)] 



(30) 



Since e <ti Q, one can average this equation over a 
timescale long compared to O -1 , but short compared to 
e _1 , yielding 

(SR 2 ){t) = £ 2 [1 + n B (uj+) + n B (w-)] cosh (et). (31) 

The fluctuations of the center-of-mass are thus expo- 
nentially amplified due to the parametric modulation of 
the magnetic field. Notice that even at zero tempera- 
ture, where initially the two bosonic modes are unoccu- 
pied [n B (uj±) = in Eq. (31)], parametric amplification 
takes place, triggered by the quantum fluctuations of the 
center-of-mass coordinate. 



IV. EFFECTS OF DAMPING AND 
ANHARMONICITIES ON THE PARAMETRIC 
RESONANCE OF THE MAGNETOPLASMONS 

In this section, we analyze the role played by dissipa- 
tion and non-linearities on the parametric amplification 
of the magnetoplasmons. As we will show, these two ef- 
fects are indeed limiting the exponential amplification of 
the bosonic modes presented in Sec. III. It is therefore 
relevant to study them as they will inevitably affect ex- 
periments. 



A. Mean-field equations of motion 

In what follows, wc adopt a different strategy than the 
one used in Sec. Ill, and focus on the time evolution of 
the averaged quantities (25) and (26) as we are primarily 
interested in the fluctuations of the center-of-mass coor- 
dinate, Eq. (29). To this end, we use the time evolution 
of the reduced density matrix p associated to the center- 
of-mass degrees of freedom, 



,P] 



\ - la 

tr=± 



(2b a pbi-bib aP -pbib a ) 



(32) 
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The first term on the right-hand side of Eq. (32) ac- 
counts for the unitary dynamics of the center-of-mass 
density matrix that evolves according to the Hamilto- 
nian (12). The second term accounts for dissipative 
processes coming from the anharmonicites [and hence 
from the coupling Hamiltonian (9)], as well as stemming 
from other sources of dissipation, such as radiative damp- 
ing, interaction with phonons, Joule heating due to the 
eddy currents generated by the electric field associated 
with the time-dependent vector potential (4), etc. In 
Eq. (32), we take dissipation phcnomcnologically into ac- 
count by assuming the Lindblad form for the reduced 
density matrix. 30 We denote 7+ and 7_ the damping 
rates for the two magnetoplasmon modes. Notice that, 
in general, 7+ 7^ 7_ due to the difference in energy of the 
two modes, as recently demonstrated in the context of 
the Landau damping of the magnetoplasmon excitations 
in metallic nanoparticles. 16 

From the master equation (32), the time evolution of 
the quantum average of an operator O is easily deter- 
mined using the identity (0)(t) = tr{p(t)0}. Intro- 
ducing the notation N = N + + N_ , AiV = N + - N_ , 
7 = (7+ + 7_)/2, and A7 = (7+ — 7_)/2, we find from 
Eqs. (12) and (32) 

N = iuj c Sf(t) (ip - ip*) - 7 7V - A 7 A7V 

\2(cp* - ip) + (b\.blblbl) + (blb\.blb + ) 

+ {b\.bt_bLb_) - (tf + b + b + b_) - (b[_b + b_b_) 

- (b + b + b_b_)] , (33a) 

-/AN ~ A-fN, (33b) 

i [u + + u- + co c 5f(t)] ip - 



AN 





i fSf(t) (N + l) 



jip - 1— (2 + 6ip + 2tp* +4N + (b\.b\.bib + ) 

+ (blblbtb_) + A{b\bH + b_) + (blblb + b + ) 
+ (bl.bl.b_b_) + 3(blb + b + b_) + 3(blb + b_b_) 
+ 2(b + b + b_b_)). (33c) 



To first order in the small anharmonicity parameter a <C 
1, the four-operator correlators appearing in Eq. (33) 
must be evaluated to order a . Since for a — 0, the 
Hamiltonian of Eq. (12) is quadratic in the bosonic oper- 
ators, one can apply Wick's theorem, yielding the mean- 
field equations of motion 

N = ioj c Sf(t) (cp-cp*)-jN- AjAN 

+ 2iaoj a [{<p - <p*) (N + 1) + (p 2 - >p* 2 ] , (34a) 
AN = - j AN - A 7 7V, (34b) 

<P = ~ i [w+ + w- + w c */(t)] <P - iy Sf(t) (N + 1) 



j<p — iaujQ {N + if + 2tp*<p + 2p 2 



+ (3(p + ip*) (N + 1) 



(34c) 



Introducing 

p(t) = exp (-i J ds [uj+ +uj- + ui c 6f(s)]j <p(t), (35) 

to first order in 77 -C 1 and within the rotating wave 
approximation close to resonance (f2 ~ uj + + we 
find from Eq. (34) 



N = e 



3 i(n-w + -w_)t ~ _j_ 



ip + c.c. 



■ jN - AjAN 



— 2aujQ 



e i(n-u+-o;_)t~ + 



<p + c.c. 



AN 





■ j AN - A-/N, 



(N+l), (36a) 
(36b) 



= ie-W-v+'^^N + 1) - 70 - 3iacj o (N + 1) 



aujQ- 



(O— CJ_|_ —UJ- )t 



(N + if + 200* 



}• 



(36c) 



B. Effect of dissipative processes on the parametric 
resonance of the magnetoplasmon excitations 

It is instructive to analyze the effect of damping alone 
on the parametric resonance of the magnetoplasmon ex- 
citations. To this end, we set the anharmonicity param- 
eter a to zero in Eq. (36). Using the initial conditions 
N(0) = n B {uj + ) + n B (w_), AN(0) = n B (ui + ) - n B (ui_), 
and 0(0) = 0, we obtain the solution at resonance 
(0 = uj + + u)_) 



N(t) 
AN(t) 



0(t) 



e" 7t 
Ve 2 + A7 2 

ee" 7 * 
v/e 2 + A7 2 

- %Kt) 

2 ^e 2 + A7 2 
+ h(t) 



h 'W e 2 + A 7 2„ 7 2^ 



A7 



{ AJV(O) 



7 e 2 + A7 2 — 7 2 



A7 



Ve 2 + A7 2 V 
7 2 - A7 2 



(AJV(O) 



2 7 e 2 + A7 2 - 7 2 



(37a) 

*?) 

(37b) 
(37c) 



where we defined the function 

2 



h{t) 



A7 2 - 7 2 



iV(0) 



sinh ( \/e 2 + A7 2 t 



7£ 



A 7 AiV(0) 



e 2 + A7 2 - 7 2 
cosh (^fe 2 + A^ 2 t 
v/e 2 + A 7 2 " 
and h'(t) denotes its derivative with respect to time. 



(38) 
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FIG. 2. (Color online) (a) Total number N(t) and (b) differ- 
ence AN(t) of bosonic excitations between the modes + and — 
as a function of time from Eqs. (37a) and (37b), respectively. 
In the figure, T = 0, A7 = -e/10, and 7 = 3e/4, \/LuTe, and 
5e/4 for the solid (red), dotted (black), and dashed (blue) 
lines, respectively. 
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FIG. 3. (Color online) (a) Total number N(t) and (b) differ- 
ence AN(t) of bosonic excitations as a function of time from 
a numerical solution of Eq. (36) at resonance (SI = uj+ + uj- ) 
for various values of the damping constant 7. In the figure, 
T = 0, A7 = — e/10, Wo = 10 3 e, Q. = i/Ecjo (corresponding to 
uj c = Wo), and a = 10 -6 . 



Thus, for e > \j "f 1 — A7 2 , the solutions (37) are ex- 
ponentially amplified, as shown by the solid (red) lines 
in Fig. 2. Indeed, from Eq. (37), we obtain for t 
(V^+Af - 7)- 1 and e > yfif - A7 2 



N(t) 



N(0) 



A7 



a/e 2 + A7 2 



+ 



AN(t) 
<p(t) 



^e 2 + A7 2 ( ^e 2 + A7 2 - 7) 
A7 



N(t), 



Ve 2 + A7 2 



2 v /e 2 +A 7 2 



A7V(0) 
(39a) 

(39b) 
(39c) 



On the contrary, for e < a/ 7 2 — A7 2 , the system reaches 
the stationary solution 

P 2 

7V st = - 



V 



AN = 



A7 2 



A 7 Ars t 



7 

(S st = — (N s 
2 7 V 



1 



(40a) 
(40b) 
(40c) 



as exemplified by the dashed (blue) lines in Fig. 2. In 
terms of bosonic excitations in the two modes + and — , 
this translates into 

o, (41) 



7t 



i.e., the detailed balance 7+-/V+ = 7 _iVI t . For the special 

case e = A/7 2 — A7 2 , we obtain from Eqs. (37) and (38), 
and for jt ^> 1 



N(t) 



AN{t) 



<p(t) 



e 

7 

A7 
27 



t - 
. 7 



27 



•JV(0) - — A7V(0) 

7 



(42a) 



e 

4 7 



r e 2 



3 

2^ 



27 



7V(0) 



A7 



AJV(0) 



(42b) 



JV(0) - — A7V(0) 

7 



(42c) 



7+ + 7- 7+7- 



i.e., a linear behavior of the solution as a function of t 
[see dotted (black) lines in Fig. 2]. 



C. Effect of anharmonicities on the parametric 
resonance of the magnetoplasmon excitations 

We now turn to the full solution of Eq. (36) , including 
anharmonic terms, i.e., a / 0. To this end, we concen- 
trate on the case where the pumping frequency is at reso- 
nance, i.e., Q = A numerical solution of Eq. (36) 
is presented in Fig. 3 for the case where e > ^/ 7 2 — A7 2 , 
i.e., when parametric amplification takes place in the ab- 
sence of the anharmonic term (a — 0), see Sec. IV B. For 
short times, both the total number of bosonic modes N(i) 
[Fig. 3(a)] and the difference between the occupation of 
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FIG. 4. (Color online) Stationary number of bosonic modes 
from Eq. (43) as a function of the driving frequency e (solid 
red line). The asymptotic behaviors for e < \/7 2 — A 7 2 
[Eq. (40a)] and e > ^/ 7 2 - A7 2 [Eq. (45)] are shown as a dot- 
ted and a dashed line, respectively. In the figure, a — 10 -6 
and W0/7 = 2 x 10 3 . 



the + and — modes, AN(t) [Fig. 3(b)] get exponentially 
amplified. This behavior is followed by a series of oscilla- 
tions, to finally reach a stationary occupation due to the 
residual (anharmonic) interaction. 

To better understand the numerical results of Fig. 3, 
we now search for the stationary solution of Eq. (36). 
Since we work to first order in the small parameters 
e/fi <C 1 and a <C 1, we neglect the terms oc ae/Q in 
Eq. (36), yielding 



1 - 



A7 : 



-j (3auj ) 2 N st (iV st + l) 2 



(e 2 + A7 2 - 7 2 ) A st - e 2 = (43) 



and 



e{N st + l)/2 
7 + 3iaw (A st + 1)' 



(44a) 
(44b) 



The solution to Eq. (43) is shown in Fig. 4 as a red 
solid line as a function of the driving frequency e. For 
e > \/ J 2 — A7 2 and a — 0, we know from our study in 
Sec. IV B that N(t) diverges exponentially [cf. Eq. (39a)], 
such that we can safely assume for that case iV st ^ 1 in 
Eq. (43). If on top of that, (e 2 + A7 2 - 7 2 )iV st > e 2 , we 
obtain 



yje 2 + A7 2 - 7 2 
3evtJo-\/l — A7 2 /7 2 



(45) 



showing that the smaller a, the more efficient is the 
parametric amplification. The asymptotic behavior of 
Eq. (45) is shown as a dashed line in Fig. 4 and com- 
pares well with the full solution of Eq. (43) (solid red 



line). For the case e < \/"f 2 — A7 2 , we find that para- 
metric amplification does not take place, and that the 
stationary number of bosonic modes is, to first order in 
a, given by Eq. (40a), as shown in Fig. 4 by a dotted 
line. 



V. EXPERIMENTAL CONSEQUENCES 

Before presenting the various possible experimental 
consequences of our proposal, we start this section by dis- 
cussing how one can trigger the parametric resonance of 
the magnetoplasmon collective modes in realistic quasi- 
two-dimensional semiconductor quantum dots. For a 
constant average electron density in the dot, the con- 
finement frequency entering the single-particle potential 
(3) is given by 18,31 



w = 



^47re e r 7 



r?N< 



1/4 



(46) 



and can be tuned by the gate voltage forming the quan- 
tum dot through the number iV e of electrons. In Eq. 
(46), eo is the vacuum permittivity, e r the relative static 
dielectric constant of the considered material, r s the 
Wigner-Seitz radius, and to* the effective electron mass. 
Using the values for a GaAs quantum dot (e r = 13, 
to* = 0.067to with to the bare electronic mass, r s = 
l.Sdg with 03 = 4irE Q e r h 2 /m*e 2 the effective Bohr ra- 
dius), and assuming N e = 10 3 , we have ujq — 1.7THz 
(huiQ = 1.1 mcV). In a static magnetic field of B = 2 T, 
the eigenfrequencies (14) are given by uj + = 5.7 THz and 
w_ = 0.5 THz. The resonance condition Q = lu + + cj_ 
for parametric amplification to occur [see Eq. (22)] thus 
requires in this case that the magnetic field is periodi- 
cally modulated at the frequency f2/27r = 0.99 THz. The 
ongoing efforts towards the production of THz sources 
should allow one to attain such pumping frequencies in 
a near future. 32-34 As we have shown in Sees. IV B and 
IV C, parametric amplification only occurs if the driving 
frequency e of Eq. (19) is larger than the linewidth 7 
of the magnetoplasmons (assumed here to be the same 
for both collective modes). With a pumping strength 
7] = 0.1, we obtain e = 0.22 THz. Assuming 7 ~ wo/10 
which is the typical value encountered in experiments, 20 
parametric amplification should occur in this case at a 
rate given by e — 7 = 50 GHz [see Eq. (39a)]. 

In what follows, we propose several ways of detecting 
the parametric amplification of the magnetoplasmons: 
First, these modes decay, among various processes, by 
radiative damping, emitting photons at the frequencies 
lj + and w_ . Triggering the parametric resonance of the 
modes with the only help of a modulated magnetic field 
should thus result in the spontaneous emission of pho- 
tons at these frequencies, without using an external light 
source, as usually employed in far infrared spectroscopy 
experiments. 

Second, as we have shown in Sec. Ill, the average fluc- 
tuations of the electronic center-of-mass coordinate are 



directly related to the to tal number of bosonic excitations 
N in the system, (5R 2 ) ~ N. As the latter is exponen- 
tially amplified up to its stationary value (45) due to the 
anharmonicity, this should result in an expansion of the 
whole electronic cloud forming the quantum dot which 
can be measured by scanning tunneling microscopy. 35 

Another measurable quantity 36 ' 37 related to the fluctu- 
ations of the center-of-mass coordinate is the magnetiza- 
tion of the quantum dot. Indeed, the magnetization op- 
erator for the electronic system described by the Hamil- 
tonian (2) reads 



M = -— Yn 

2m* ^ 

i—\ 



[Pi + eA(rj,t)] 



(47) 



Using Eq. (4) and introducing center-of-mass and rela- 
tive coordinates as in Sec. II, the magnetization operator 
separates into M = M cm + M ro i, with 



2m* 



N c eB(t) R2 



(48) 



the contribution from the center-of-mass dynamics, and 
M rc ^-^-Efc + ^ 2 )e, (49) 



the one from the relative coordinates. Since the latter is 
induced by a bath of fcrmionic quasiparticles, the mod- 
ulation of the magnetic field does not lead to their para- 
metric amplification. In fact, the only result of the mod- 
ulation is a small periodic variation of the relative mag- 
netization that averages to zero due to the fast modula- 
tion. As a consequence, the average magnetization result- 
ing from the relative coordinates is constant in time and 
equals its value without parametric modulation, (M rc i). 

In contrast, the magnetization (48) resulting from the 
center-of-mass dynamics can be conveniently rewritten in 
terms of the bosonic operators (10). Taking its expecta- 
tion value, and averaging over fastly oscillating terms, we 
find, to first order in the parametric modulation strength 
V, 



(M cm )(t) = — /ige z < AN(t) 



LOc/2 



[N(t) + 1] 



(50) 

with /ig = eh /2m* the effective Bohr magneton (fi^ = 
0.86meV/T for GaAs). As both the total number of 
bosonic excitations N(t) and the difference AN(t) are 
exponentially amplified due to the parametric modula- 
tion to finally reach stationary values [see Eqs. (44a) and 
(45)], the resulting magnetization of the dot should dra- 
matically increase when the parametric amplification of 
the magnetoplasmons takes place as compared to its equi- 
librium value. 



VI. CONCLUSION 

In conclusion, we have studied the possibility of para- 
metrically amplifying bosonic collective modes in finite- 
size fermionic systems. Specifically, we have shown that 
the magnetoplasmons in quasi-two-dimensional semicon- 
ductor quantum dots can be parametrically amplified 
by modulating the magnetic field perpendicular to the 
nanostructure. Moreover, we have demonstrated that 
damping mechanisms and anharmonicities of the confine- 
ment lead to a saturation of the parametric resonance. 
We have further discussed the implementation of our 
proposal in realistic experimental samples and suggested 
measurements that should present clear signatures of the 
parametric amplification of the magnetoplasmons. 

Our predictions could in principle also apply to the 
magnet oplasmon modes in metallic nanoparticles. How- 
ever, as the damping of these modes is much stronger 
than in quantum dots, 16 it may be much more difficult to 
trigger the parametric amplification in metallic nanopar- 
ticles. 
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Appendix A: Parametric resonance of a single 
harmonic oscillator: classical vs. quantum treatment 

In this Appendix, we consider a single harmonic oscil- 
lator of mass m whose frequency ui(t) = uiq[1 + r\ sin (fit)] 
is periodically modulated in time at the frequency f2. 
Denoting q and p the position and the momentum of the 
oscillator, its Hamiltonian reads 



H osc (t) = 



P 

2m 



(Al) 



In what follows, we assume that the strength rj of the 
periodic modulation is much smaller than one, allowing 
for a perturbative treatment of the Hamiltonian (Al). 
For an arbitrary r], the Floquet formalism would be 
the appropriate one. 38 Before we tackle the quantum- 
mechanical treatment of the Hamiltonian (Al), we recall 
in the next section the classical case. 



1. Classical treatment 

In order to find the classical trajectories associated to 
the Hamiltonian (Al), as well as to stress the analogy 



9 



with the quantum mechanical treatment of the paramet- 
ric resonance presented in the Sec. A 2, it is convenient 
to introduce the new variables 



z = (yjmojoq + 1P 
V2i V y/mwo 

- 1 ( I [p 

z = —= y/mu} q , 

V2i V \/mu} 



(A2a) 
(A2b) 



Since z and z are canonically conjugated, i.e., the Poisson 
bracket {z,z} — 1, the Hamilton equations of motion 
read z = —dH osc /dz and z = dH osc /dz. With Eq. (A2), 
the Hamiltonian (Al) transforms into 

H OBC (t) — iw [1 + 77 sm i^t)] zz + vq — sin (Q£) (z 2 + z 2 ) , 

(A3) 

yielding 

z = —iluo [1 + i] sin (fit)] z — iujoV sin (Qt)z. (A4) 
Introducing 



2. Quantum treatment 

At the quantum level, the Hamiltonian (Al) is conve- 
niently rewritten as 

H osc (t) = huj [1 + 77 sin (fit)] (tfb + 

+ 7 ? / ^sin(r!i)(V 2 + ^ (All) 

in terms of the lowering and raising bosonic operators 
1 ( q , ipt, 



b = 
&t = 



V2 
1 

T2 



u osc 



\p£ Q 



(Al2a) 
(A12b) 



respectively. Here, £ osc = ^h/mojQ denotes the oscilla- 
tor length. The Heisenberg equation of motion for the 
operator b reads 



z{t) = exp ^-io; J As [I + r)sm{ns)}^jZ(t), (A5) b = -iu [I + rj sin (fit)] b - iu) r) sin (flt)tf 



(A13) 



Eq. (A4) becomes, to first order in 77, 

Z = -iu r)sin(flt)e 2iUot Z. 



(A6) 



and is similar to the classical equation of motion (A4) for 
the classical variable z. It can thus be solved in exactly 
the same way, yielding 



Close to the resonance condition f2 ~ 2wo, we keep only 
the dominant term cx e ~H n ~ 2w °)* such that Eq. (A6) de- 
couples and yields 



b(t) 



-iujQt 



n. 



z + i(n-2Lu Q )z + {^f) z = o. 



(A7) 



i^o (e ia +t 



+ e in -' 



)&(0) 



fot(O) 



(A14) 



With the initial conditions Z(0) = z(0) and Z(0) — 
7]ujqZ(0)/2, we finally obtain 



z(t) = 



-ILUot 



with 



2 V 



ft_e if2+t - o+e 



*(0) 



**(0) 



(A8) 



n± = w - 'j ± ^ - 2w ) 2 - (^o) 2 . (A9) 

Thus, when |fi — 2wo| < rjLUo, the frequencies f2± become 
imaginary, which leads to the exponential amplification 
of the motion of the oscillator. Indeed, coming back to 
the original variables of the Hamiltonian (Al) and as- 
suming f2 = 2loq, one obtains 

q(t) = q(0) cos {uo a t)e^ Qt / 2 + A sin ( Wo i) e -"^*/ 2 , 

itlujq 

(AlOa) 

p(t) = p(Q) cos (w i)e-" w °*/ 2 - mcj g(0) sin {Lj t)e nuj ° t/2 . 

(AlOb) 

Notice that the parametric amplification only takes place 
within the above approximate treatment if initially, 
q(0) ? 0. 



with the frequencies f2± given in Eq. (A9). Thus, at the 
quantum-mechanical level, the lowering and raising oper- 
ators are exponentially amplified whenever the resonance 
condition |fi — 2cjo| < t/wo is fulfilled. 

Considering that the harmonic oscillator is in its 
ground state before the parametric amplification takes 
place (i < 0), and assuming ft — 2u; , the average num- 
ber of bosonic excitations 

(6t6)( t ) = I[coBh(» W ,t)-l] (A15) 

exponentially increases for increasing time. During that 
process, the average position and momentum arc both 
vanishing, (q)(t) — (p)(t) — 0, while the corresponding 
fluctuations are exponentially amplified, 

i 2 

(Sq 2 ) (t) = -^p [cosh (r]uiot) + cos (2u> t) sinh (rjLJot)] , 

(Al6a) 

[cosh (rjUJot) — cos (2uj t) sinh (rjwot)] . 



(Sp 2 )(t) 



2P 



(Al6b) 



The quantum fluctuations of the oscillator in its ground 
state thus provide a seed for the parametric resonance in 
the quantum case. 
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Appendix B: Parametric resonance of the 
magnetoplasmons: classical treatment 

Classically, the center of mass Hamiltonian (7) in the 
absence of the quartic correction (A4 = 0) can be solved 
by a series of two successive canonical transformations. 
First, we introduce 



£=-L(X + iY), P e = ±(P x -i Y ), (Bla) 
C = ^=(X-iY), P e , = -L( Px + iY), (Bib) 

with Poisson brackets {£, Pj} = {f*,P£*} — 1: such that 
the Hamiltonian (7) transforms into 



H an (t) 



M 



+ M Wq + 



(B2) 



Second, we introduce the canonically-conjugated vari- 
ables 



z+ 



V2i 


he- 




V2i 







/3 



(B3a) 



(B3b) 

with /3 = -\/^( w o + w c/4) 1/ ' 2 an d the Poisson brackets 
{z±,z±} = 1. In terms of these variables, the Hamilto- 
nian (B2) reads 



+ i^6f(t)(z + z-+z + z-), (B4) 



where the function Sf(t) is defined in Eq. (13) and the 
frequencies uj± are given in Eq. (14). The Hamilton equa- 
tions of motion associated to the Hamiltonian (B4) thus 
read 

z± = = -iw± [1 ± Sf(t)} z ± -i^Sf(t)z T . (B5) 

Being similar to the quantum mechanical Heisenberg 
equations of motion (15) for the bosonic operators b±, 
the above equation can be solved in exactly the same 
way as in Sec. III. Going back to the original variables of 
the center of mass Hamiltonian (7), we find that close to 
resonance [\Sl — u + — <C e with e given in Eq. (19)], 
the classical motion of the center of mass is described by 



the trajectories 

X(t) = {X(0) [cos (u+t) + cos (w-t)] 

-Y(0) [sin (w+t) - sin (u)-t)]} 
+ {Pxr(0) [sin(w + <) +sin 



e et/2 

~2~ 



+Py(0) [cos (w+t) - cos (w_t)]} — 5-, (B6a) 



e -et/2 
~2/3 2 



3 et/2 



= {F(Q) [cos (cj+i) + cos (oj-t)] 

+X(0) [sin - sin (w_f)]} 

+ {Py(0) [sin + sin 

e -ei/2 

-Px(0) [cos ( w+ t) - cos (w_t)]} -—a". (B6b) 

In Fig. 5, we show typical trajectories of the center 
of mass from Eq. (B6) in absence of parametric modu- 
lation [Fig. 5(a)] and with parametric modulation [Fig. 
5(b)]. Without parametric modulation, the center of 
mass has trajectories corresponding to a two-dimensional 
harmonic oscillator in a perpendicular magnetic field 
[Fig. 5(a)]. The corresponding radial coordinate shown 
in Fig. 5(c) thus regularly oscillates between and its 
maximal amplitude. In contrast, in the presence of para- 
metric modulation [Figs. 5(b) and 5(d)], the motion of 
the center of mass is exponentially amplified. 



Appendix C: Toy model: Parametric modulation of 
two coupled fermionic modes 

In this Appendix, we consider two coupled fermionic 
modes with frequencies e + and e_ whose corresponding 
Hamiltonian is parametrically modulated. In the sequel, 
we show that due to the fermionic nature of the modes, 
Pauli principle prevents the parametric amplification to 
occur. 

As a toy model, we consider a Hamiltonian similar 
to Eq. (12) without quartic correction, except that the 
bosonic operators in Eq. (12) are replaced by fermionic 
ones, 

Htoy(t) = ^ he a [1 + a?/ sin (Sit)] S a c a 

+ 77^ sin (Sit) (c{cl + c_c+) . (CI) 

Here, c\ / c a creates/annihilates a fermion in mode = ±. 
In Eq. (CI), 77 <C 1 is the strength of the parametric 
modulation at frequency SI and A is the frequency scale of 
a broken symmetry term. The corresponding Heisenberg 
equations of motion read 



-ie± [1 ± 77 sin (Sit)] c± =p ir?— sin (Orj)c^. (C2) 



A 

~2 
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FIG. 5. (Color online) Upper panel: Typical trajectory (X(t),Y(t)) of the center of mass coordinate (a) in the absence of 
parametric modulation (77 = 0) and (b) with parametric modulation (77 = 0.1) for times up to £ ma x = 250/o;o. Lower panel: 
radial coordinate R(t) = \JX 2 {t) + Y 2 (t) as a function of time for (c) 77 = and (d) rj — 0.1. For all trajectories, the initial 
conditions are X(0) = Xq, Y(0) — and Px(0) = Py(0) ~ 0. In the figure, the cyclotron frequency u) c = ujo and the pumping 
frequency is at resonance, i.e., = oj+ + 



Along the same lines as in Sec. Ill, introducing the new 
operators 

c± (t) = exp ^-ie± J ds [1 ± rjsin (£ls)]j c±(t), (C3) 

we obtain for Eq. (C2) to first order in rj and within 
the rotating wave approximation, which is valid for Q ~ 



with 



c± + i (ft - e + - e-) c± 



c± = 0. 



(C4) 



Solving for the latter equation leads to 

e -ie±t/fi 



c±(t) = 



±i 



n_e in+t - n + e lU ~ l ) c±(0) 



Q-t\ 



!?^(e^+t-e in -*)4(0) 



(C5) 



e_ - Q 1 

± - 

2 2 



(C6) 

Thus, contrarily to the bosonic case [cf. Eq. (22)], the 
frequencies fi± are real for any pumping frequency f2, 
such that the fermionic operators (C5) are not paramet- 
rically amplified. This is a direct consequence of Pauli's 
principle. 
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